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Introduction 


A proton of low energy moving in the equatorial plane of the earth 
will experience drift motions due to both the magnetic field (magnetic 
gradient drift only, if the field is assumed to he that of a dipole) and 
the electric field. The electric drift again separates into two parts - 
the drift due to the "main" electric field (or 1T convection electric field 11 ) 
existing in the frame of the earth, and that due to the earth T s rotation. 


Such motion has heen analyzed hy Chen \ 1970 1 who assumed that the 
main field in the equatorial plane is a constant field from dawn to dusk. 

He obtained a somewhat complex variety of proton orbits: the main purpose 
of this work is to classify such orbits and to devise a simple procedure 
by which one can determine, in the energy spectrum observed by a satellite 
in the equator, which energies (of equatorial particles) correspond to 
trapped trajectories and which to open ones. The elctric field used here 
f Stern 3 1974 1 differs somewhat from the model used by Chen, giving a 
more realistic fit to high-latitude observations of the electric field, 
but all results obtained here could be easily generalized for Chen r s model. 


One result indicated by this work is that at distances of 4 - 6 
earth radii, a transition from trapped proton orbits to open trajectories 
leading to the tail occurs at about 10 kev, the precise value depending 
upon local time. Such a transition also seems to be indicated by particle 
observations of Smith and Hoffman \ 1975a » bj , using Explorer 45* The 
energy spectrum (at magnetically quiet times) of equatorial protons above this 
energy can be explained by charge exchange [ Smith and Hoffman , 1974 j , but 
increased flux observed below it seems to be related to the influx of parti- 
cles on open orbits from the tail. 


Qualitatively the situation may be viewed as follows. At very low ener- 
gies and near enough to earth, electric drift due to the co-rotation field 
is dominant, and charged particles tend to co- rotate and are thus trapped. 
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At. relatively high energies all electric drifts may be neglected and guiding 
center motion is mainly due "to : » i c drifts "which also tends to keep 
^joj^icles trapped in the vicinity of the earth. 

Nov for equatorial electrons these two drifts are in the same direction 
and tend to reinforce each other, while for equatorial protons the directions 
are opposed o There may thus exist for protons an intermediate energy range 
in which these two drifts approximately cancel and the drift in the main electric 
field ( "convection field”) dominates the motion: drift orbits due to this field 
are open and generally follow the no on* - midnight direction, allowing an expla- 
nation for the previously mentioned features oi the energy spectrunio 

The Electric Field 


1 that the 

Following previous work T Stern , 197^ J Ve assume ^ main electric field 
in the rotating frame has a form derivable from a scalar potential* On 
closed field lines our model gives 


E = - 

k/2 

<j) = - (j)^ ( c4 q / o£ ) sin Y 


(1) 

( 2 ) 


Here ol, is an Euler potential for the dipole magnetic field (which is the 
magnetic configuration assumed), given by 

oC = a g° $ ^ sin^ & (3) 

1 

where g° = - 3.1 * icf^ MKS is the dipole term in the harmonic expansion 
of the earth's magnetic field, a is the earth’s radius and f - r/a . 
Also, Q is the value of oL on the boundary of the polar cap and 2 <j> Q 
is assumed to be the voltage across the cap. The power k reflects the 
rate with which j E | decreases near the polar cap boundary: Chen's model 
is obtained if k = 2 , but electric field measurements at high latitudes 

j^ Heppner , I 972 J suggest that k = 4 is more realistic, and this value will 
be assumed here . 
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The total electric field - with axisymmetric rotation added - is 
E = - V(^ + M 

with $ = a (h c4 ( 5 ) 

2 

where U) = 7„3 • 10~^ rad/sec is the earth's angular velocity. Thus 
in the equatorial plane , with k = 4 

<j> = - (h ( oC /a. g°) $ sinY + a 2 g° Cd / ? (6) 

/ IQ O 1 X 

The Hamiltonian 


The Hamiltonian for the motion of a particle of charge e in the dipole 
equator, in the presence of the fields described here, is 

H = p B + e <f> ( 7 ) 

In dealing with protons, if e is dropped and 4 is reckoned in volts, then 
all energies are in electron volts and the magnetic moment p of the par- 
ticle (which is assumed to he conserved) is conveniently expressed in ev/B . 
The magnetic field B itself is given by 

B = |g°| /? 5 (8) 

The Hamiltonian is a constant of the motion - that is, contours of 
constant H in the ( 9 > ¥ ) plane for any given describe trajectories 
of protons having that value of p * 

For convenience we may use the constancy of H T = H/jg°| > which can 
he expressed in terms of positive constants q and p (no relation to 
canonical variables) 



" 5 " 


H* = H - q. J 2 sinT - p$ 1 

(9) 


v here 

* = " i ( ^o /a /( 6 °| 

(10) 


P = a 2 u) = 2.9 ■ 10^ m 2 /sec 

(11) 


If the radius of the polar cap is taken as a/t and the 

mean electric 

field there is E volts/m , then 

0 



q = 2 * 10 8 E 

* O 

(12) 


The calculation may he made more general if as many constants as possible 

are removed by suitable normalization of quantities* To remove q 

and p 

one first divides ( 9 ) by q and then defines new units 
that ? is replaced by 

of length. 

such 

R ■= ? (q/p) 1 / 5 

(13) 


Specifically, if one defines 



W = H '/( q p 2 ) 1 / 5 

(Ik) 


M = p (q 2 /? 5 ) 1 ^ 5 

(15) 


then equation ( 9 ) transforms into 



¥ = ME' 5 - B 2 sin N* - R _1 

(16) 



which will he the basic conservation law used in what follows. 
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Qualitative Analysis 


If W is considered as a function of (r, 'f ), proton trajectories in 
the equatorial plane will be given by contours of constant W and "last 
closed trajectories" will be characterized by contours which cross them- 
selves (see Figures 1, 2 and 3) • At a crossing point the direction of V W 
is not defined and therefore 


yw = 0 (17) 

The reverse however is not true : V W may vanish at a point which is 
not a crossing point but rather the limit point of a set of nested closed 
trajectories (e.g. point in Figures 1 and 2). 

• The -derivative of (16) shows that (17) may only be satisfied on the 
dawn-dusk line, at points where 

sin 'f = +1 (18) 

It is instructive to examine the dawn side and the dusk side of that line 
separately. 

On the dusk side (sin^f = l) the vanishing of the R-derivative of 
(16) gives 

1 = 3 M R -2 + 2 R 5 = |(R) (19) 

The function ^(r) tends to infinity at both 0 and 00 and it 

has a single minimum at R = o The value of £{r) at the minimum 

point will he smaller or larger than 1 > depending on whether M is less 
or more than 

M q = = 0.068if- (20) 

and accordingly (19) - and hence (17) - will be satisfied either at two 
points or at no points at all. 
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On the dawn side ( 17 ) requires 

1 = 3 M R " 2 - 2 R 5 = *\,(R) (21) 

The function v»(r) has a negative derivative and descends continuously 
from ^( 0 ) = 00 to = - ©0 . There exists therefore on the 

dawn side a single solution of ( 17 ) for a ll values of M • 

The values of W corresponding to solutions of ( 17 ) will he called 
critical values . In particular, Wi and W 2 will denote critical values 
corresponding to the dusk side (Wj having a larger value of R than W 2 ) and 
W 3 will denote critical values on the dawn side. In all cases the critical 
values depend on M , so that W a , W 2 and W 3 are all functions of M , displayed 
in graphical form in Figure 4 . 

Interpretation 


If the point for which W was derived is a crossing point, then the 
contour W(M, R, 'f ) = W ± (m) , for any given M, gives in part the boundary 
between a region of trapped trajectories and a region of open ones. If, 
on the other hand, represents a limit point of a nested set of trapped 

trajectories, then the line W(M,R, Y ) = W (m) shrinks to a point and does 
not properly exist. It will turn out that the first case applies to and 

W 3 while the second holds for W 2 , 

As M is varied, the map of contours of constant W(M, R, 'p ) gradually 
changes. Three regimes then exist and will be denoted by numbers 1, 2 and 3 > 
they are illustrated by the Figures bearing the same numbers and by the 
numbered ranges of M in Figure 4 , In regime 1 (illustrated schematically 
by Figure 1 and also, by Figures 2 in the work of Chen f 1970 ] ) w l > w 3 

and the trapped region defined by encloses the one vhich is defined by 

W 3 , so that it alone determines vhether an orbit is trapped or not. In 
regime 2 (Figure 2 ; also Figures 3b and 4a in the "work of Chen [±2Z2] > 

where Figure 3a gives the transition between regimes 1 and 2 ) < W 3 and 

the trapped regions determined by Wj. and W 3 are disjoint. Finally, in 
regime 3 (Figure 3 here, and Figures 4b and 5 Chen‘s vork) the trapped 
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region defined by Wi disappears, since ( 19 ) no longer has any real solutions, 
leaving just the trapped region defined by W 3 , which in all three regions 
surrounds the origin. 

It is instructive to sketch qualitatively the variation of W along 

the dawn-dusk line, and this is done in the lower parts of Figures 1, 2 

and 3 • Note that this variation has extrema at all points at which ( 17 ) 

are 

holds, although the shapes of these extrema only schematically shown 
in the figures. 

Each of the contours corresponding to a crossing point in one of these 
figures - except for the W 3 contour of Figure 1 - consists of two parts 

- one of them forming a closed loop, the other extending to infinity at 

both its ends. These will be called the closed and open branches of the 
contours of or W 3 . 

It is easily shown that the crossing point - which belongs to both branches 

- is the point closest to the origin on the open branch and the point farthest 
from the origin on the closed one. The proof is as follows. For given 

fixed values of and M , these contours satisfy the equation 

W j _ = M R" 3 + R 2 sin'f - R x (22') 

(plus for i = 3 , minus for i = l). They therefore define a functional 
relationship R = R( v f ) 9 with an inverse function ¥ = ^(R) ♦ At the point 
on the contour at which R is extremal dR/d^ vanishes and therefore 
d'f/d R diverges. From (22), however, 

± cos M* ay/dR - 2 W i R" 3 - 5 M r" 6 + 3 R~ 4 ( 23 ) 

Therefore such extremal points only occur in the dawn - dusk plane 
where cos Y vanishes. There exists only one such point on open branches, 
namely the crossing point, which therefore is a point of minimum R . 
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In the immediate vicinity of the crossing point the contours of 
constant assume the form of an X , and the curve of constant R 

through that point is a circular arc tangential to a line which bisects 
that X • One side of the X then belongs to the open branch and R is 
minimal there, from which it follows that for the other side, belonging to 
the closed branch, R is maximal at the crossing - point. However, the 
closed branches (except for that of W 3 in Figure 1 ) have only two 
points with cos 'f = 0 . The other point with that property is therefore 
a minimum point for R , showing that R at the crossing point is an 
absolute maximum value for the entire closed branch® 

The Profile 'Line 


Suppose one is observing equatorial protons with various energies at 
some given point (r, ‘f ) at which the magnetic field intensity is B . 
Given the various constants of the model one can then derive R for the 
observation point and values of M for each observed value of the kinetic 
energy K • If S '•'he distance in earth radii and p is the magnetic 

moment in kev /gamma, then 


0.412 E 1 ' 3 <? 

0 ^ 

(24) 

2 /3 

58.4 E q ! u 

(25) 


The objective of this section is to develop in that case a simple way 
of distinguishing which values of M (or of K) correspond to open 
trajectories and which to closed ones trapped near the earth. 

For every value, of observed energy there corresponds by (25) a value 
of M , and by equation (l6) the corresponding value of W can then 
be computed. Since (R, v f ) are fixed, equation (16) shows that all related 
pairs of values (W, m) observed at the given point belong to a single 
straight line in the (W, M) plane , which will be called the profile line 
of the observation point. The type of trajectory may change from trapped to 
free or vice versa only at points at which the profile line crosses one of 



10 - 


the lines 

W = W X (M) , W = w 3 (m) (26) 

To establish which crossings do represent a change in type of trajectory 
we follow the profile line from its right-hand side, corresponding to large 
values of M • Equation (21 ) shows that for any given M the normalized 
radial distance R 3 to the crossing point of W 3 satisfies 


3 M = R 3 2 + 2 R 3 5 

In the limit of large 


M = 2 R 3 5 /3 
W 3 - 5 R 3 2 /3 


(27) 


(28) 


On the dusk side the line W = W 3 then approaches the origin "within a 
distance R , vhere £ ^ 1 (see Figure 3) • Substituting in (l 6 ) 

gives, in the same limit 

2 - 3 S 5 = 5 $ 3 

4 5 5 2 * 5 (29) 

The inequality yields 

5 > 3" ° A (30) 

meaning that the normalized minimum distance R on the closed branch of 
W = W 3 also becomes large under these circumstances. Since the value of 
R for the observation point is fixed and finite, it follows that the 
upper end (M of the profile line alway s starts at points correspon- 

ding to the interior of the closed branch in Figure 3 • 
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At M = 0 equation (27) shows that R3 0 , i.e. the observation point 
+.*,« ,-™ 0 r ninsprl branch of W-a . Thus the profile line crosses 

j_d uuuoxllc ^ — — w - kind 

the line W = W 3 (m) at least once. If the first crossing of this occurs 
in regimes 2 or 3 , it changes the type of trajectory from trapped to open. 
If the change occurs in regime 1 , no change occurs, because all crossings 
of W3 there are from one trapped mode to another (see Figure l) . 

A crossing of the line W = Wj_(m) will change the status of the 
*kra jectories only if the closed branch of that line is involved. This 
will generally be the case, but a simple check is at any rate available to 
ascertain this, based on the theorem proved in the preceding section. To 
perform this check, one notes the value of M associated with the crossing 
point and then compares the value 01 R associated with this value of M 
by equation (21 ) (Table 1 ^Figure 5) with the value of R derived for 
the observation point. If the latter is smaller, the crossing is on the 
closed branch. 


It is interesting to examine the family of profile lines having the 
same value R Q of the normalized radial distance but different values of 
magnetic longitude 'f . All these lines are parallel and are located 
inside the strip bounded by 

W = MR 0 ‘ 3 - R 0 2 + 

and 

W = M R 0 -3 + R 0 2 + 

For R = R q there will exist some positive solution for M satis- 

fying equation (21). The point on the profile line (32) corresponding to 
Mq then also belongs to the line W = W 3 (m) > since both conditions (l6) 
and (21) are fulfilled. In fact, the profile line will be tangential to 
the W 3 curve, since (it can be shown) both lines have the same slope at 
Mq . Similarly, the bounding line (31) is tangential to W = Ws(m) , pro- 
vided the two lines meet. The situation is illustrated in Figure 6 , for 
a radial distance of 5 earth radii and given E Q - 0.01 volt/m • A third 
line in the middle of the strip represents observations on the noon-midnight 
meridian, and heavily drawn portions of the profile lines represent open 
trajectories . 


R (31) 

o 

Rq" 1 (32) 



12 


Conclus ions 


At distances of the order of 5 earth radii, there will exist for each 
observation point a range of proton energies for which trajectories are 
connected to the tail, while at lower and higher energies trajectories 
are trapped* 

Smith and Hoffman [ 1975a , Fig. 5 ] observe at quiet times an energy 
spectrum of equatorial protons which has two peaks, separated by a dip in 
the range 10 - 100 kev. The dip and the high energy peak can be explained 
Smith and Hoffman 1975b , 1974 ] in terms of charge exchange effects, but 
the peak at the lower energies falls in the region in whic^present theory 
indicates open trajectories and therefore appears to represent particles 
entering from the tail; the source being intense enough to overcome 
charge-exchange losses o 

The range of open trajectories is also expected to have a lov~energy 
boundary and this can be compared to a similar boundary observed by 
Mcllwain [ 1972 , Figure 6 ] , shifting towards lower energies as local 
time varies from evening to early morning hours. 

This shift agrees with the predictions of the model, but there exist 
notable differences the most significant being the absence of a similar 
gff 0 ct, on the day side. However, plasma sheet particles require relatively 
long times to reach the day side, during which their motion may be modified, 
and in addition, a better model of the magnetic field is probably needed 
at 6.6 earth radii, where Mcllwain's observations were made. 

The actual electric field of the earth varies in time and therefore, 
strictly speaking, it cannot be represented by a scalar potential. If however 
its variation is slow and consists mainly of a change in E q , then to lowest 
approximation the method may still be used - at least on the night side - and 
only the value of used in (24) and ( 25 ) is changed. Note that the 

variation of R derived from (24) implies a time-dependent profile line, 
which might cause the trapping of particles initially on open trajectories. 
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Two final points may be worth pointing out. First, it should be noted 
again that the basic reason protons have more complex trajectories and can 
reach points near earth more easily than electrons of the same magnetic 
moment is that their co- rotation drift opposes their magnetic drift. If 
the earth’s magnetic axis were reversed, these properties would hold for 
electrons rather than for protons. 

Now from Pioneer 10 observations it appears that the magnetic dipole 
axis of Jupiter indeed is. aligned in the opposite way from that of the 
earth, relative to the rotation vector. Observations made by that space- 
craft also indicate that Jupiter's magnetic field contains a very sizable 
ring current - causing distortions of the field comparable to those found 
in the earth's magnetic tail - and that in the trapped radiation round 
there, electrons greatly predominate (by contrast, the terrestrial ring 
current is mainly due to protons). This suggests that Jupiter might have 
an electric field similar to the earth's and that the same processes 
described here for protons are responsible for Jupiter's trapped electrons. 

The second point is an analogy between the method of the profile line 
and the Clemmow-Mullaly -Allis (CMA) diagram for the propagation of waves in 
a cold two-component plasma ^Stix, 1962 ^ . In that case, the propagation 

of waves in any given plasma depends on an appropriate dispersion relation 
and there exists an infinity of such relations., depending on the magnetic 
field intensity B and on the plasma density n o If, however, one merely 
asks 'which propagation modes are possible in a given plasma, one only needs 
to knov under vhat conditions modes appear or disappear Q In the plane in 
which the cartesian coordinates are B and n (or suitable functions thereof) 
these conditions generally occur along certain lines, and a glance at the 
diagram in which all such lines are drawn, dividing the plane among various 
modes - which is the CMA diagram - quickily yields a large amount of useful 
qualitative informat i on „ In this analogy. Figure (4-) corresponds to the 
CMA diagram, while contoured planes like those shown in Figures 1 - 3 , 
in which lines of constant W are mapped, correspond to graphs of the 
dispersion relation* 
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tions to Figures and Table 


Schematic structure of contours of constant W in 
the equatorial plane, for values of M for which W x > W 3 
( M < 0.04042 ) . Below the contours the variation of W 
along the dawn-dusk line is schematically shown. 

Similar to Figure 1 , for the range 0.04042 > M> M q . 
for which W 3 

Similar to Figure 1, for M > M q , the regime in which 
the W x contour no longer exists. 

The contours of W x (m), W 2 (m) and Ws(m) in the (W,M) plane. 

The relation between M and R along the line W = W x (m) . 

Similar to Figure 4 , with the profile lines drawn for 
? = 5 , sin "f = 0 , + 1 . 


Values of M corresponding to various values of R along 
the line V = W x (m) . The plot of this relation is given in 
Figure 5 • 


Table 1 





Table 1 


r M x 100 


R 


0.585 

6.840 

0.6 

6.816 

0.62 

6.706 

0.64 

6.495 

0.66 

6.171 

0.68 

5.720 

0.7 

5.129 

0.72 

4.381 


0.73 

0.74 

0.75 

0.76 

0-77 

0.78 


M x 100 

3.943 

3.it60 

2.930 

2.350 

1.718 

1.032 

0.290 


0.79 
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